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Abstract: Quantum mechanical models with extended supersymmetry find interesting 
applications in worldline approaches to relativistic field theories. In this paper we con- 
sider one-dimensional nonlinear sigma models with O(A^) extended supersymmetry on the 
worldline, which are used in the study of higher spin fields on curved backgrounds. We cal- 
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expansion, using both canonical methods and path integrals. The latter are constructed us- 
ing three different regularization schemes, and the corresponding counterterms that ensure 
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1 Introduction 

Quantum mechanical models with extended supersymmetry find useful applications in 
the worldline description of relativistic field theories. Indeed fields with spin S can be 
described in four dimensions by quantizing particle actions with N = 2S extended local 
supersymmetry on the worldline [1-3]. While complete in four dimensions, these models 
describe only conformal invariant fields in other dimensions [4, 5]. They can be consistenly 
defined not only in Minkowski space, but also in maximally symmetric spacetimes [6] and, 
more generally, in conformally flat spacetimes [7] . Upon gauge fixing, the worldline actions 
of the spinning particles moving on such spaces give rise to an interesting class of one- 
dimensional nonlinear sigma models possessing extended rigid super symmetries. The goal 
of this paper is to analyze the corresponding quantum mechanics. 

In particular, we are interested in computing the transition amplitude for the O(A^) ex- 
tended supersymmetric nonlinear sigma models in a euclidean short time expansion since, 
in applications to higher spin field theory, the transition amplitude can be used to study 
ultraviolet properties of propagators and one- loop effective actions [8]. We will achieve this 
result using two different methods. The first method employs canonical quantization, and 
starting from the operatorial definition of the transition amplitude we compute it perturba- 
tively in the euclidean time /3 by using the commutation properties of the various operators. 
The final result yields a perturbative solution of the heat equation (the Schrodinger equa- 
tion with imaginary time) and identifies a benchmark for equivalent calculations in terms of 
path integrals. This canonical approach has been already employed in [9] for the = 0, 1, 2 
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supersymmetric quantum mechanics, that can be defined on any curved manifold (see [10] 
for a review of the method in the bosonic case). We extend that computation to arbi- 
trary N for the O(A^) extended supersymmetric nonhnear sigma models. For N > 2 the 
extended supersymmetry may be broken by a generic curved target space, though it can 
be maintained on locally symmetric spaces. Nevertheless, we use an arbitrary target space 
geometry since for the present purposes we do not need to gauge supersymmetry, and can 
view the latter just as an accidental symmetry emerging on particular backgrounds. 

The second method we employ for computing the transition amplitude makes use of 
path integrals. Our reason for considering this approach is that, in typical first-quantized 
applications, canonical methods are used first to identify the precise operators of interest 
and path integrals are used next to perform more extensive calculations and manipula- 
tions. A classical example is the calculation of chiral anomalies and the proof of index 
theorems by worldline methods [11-13]. It is therefore important to be able to reproduce 
the transition amplitude with path integrals. Path integrals for particles moving on curved 
spaces can be quite subtle, and their consistent definition needs the specification of a reg- 
ularization scheme which carry finite counterterms to guarantee scheme independence, see 
[10] for an extensive treatment and [14] for a short discussion. We will use three differ- 
ent schemes, for completeness and because each one carries its own advantages. The first 
scheme, time slicing (TS), can be deduced directly by using operatorial methods [15, 16] 
and can be related to the lattice approximation of the propagation time. A second scheme, 
mode regularization (MR), is related to a momentum cut-off (or, more properly, energy 
cut-off in quantum mechanics) and allows the introduction of a regulated functional space 
to integrate over [17-21]. Finally, dimensional regularization (DR) is defined as a purely 
perturbative regularization but has the advantage of carrying only covariant counterterms 
[22-26]. For each of the three regularization schemes we find the corresponding coun- 
terterms that ensure scheme independence, making them ready for extending to curved 
backgrounds the worldline approach to higher spin fields initiated in [27]. 

Let us now describe the precise form of the (supersymmetric) quantum mechanics we 
are interested in. We consider a particle moving in a curved space Ai of dimensions D and 
metric ^^^y. It is described in phase space by bosonic coordinates and momenta {x^^p^)^ 
where /x = 1, . . . , is a curved spacetime index, and by Majorana worldline fermions (i.e. 
real Grassmann variables) '0^, where z = 1, . . . , is a flavour index and a = 1, . . . , is a 
flat spacetime index related to curved indices by the vielbein e^, defined by g^y — rfab^^^t 
with r]ab the flat tangent space metric (which can be taken either minkowskian or euclidean, 
according to the desired applications). Quantum mechanically the bosonic variables satisfy 
the usual commutation relations [x^,piy] = iSH^, and the fermionic ones give rise to a 
multi Clifford algebra {j/j^^j/jj} = rf^^Sij. The general hamiltonian operator we wish to 
consider involves a free covariant kinetic term Hq, a four- Fermi interaction depending 
on the Riemann tensor and carrying a coupling constant a, and a contribution from an 
arbitrary scalar potential V which depends only on the spacetime coordinates (and 
which in most applications is proportional to the curvature scalar R) 

H^Ho + aRabcdrA^j^'^j + V (1.1) 
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where 
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-UJ^abi^fi^i 



(1.2) 



with uji^ab the spin connection and ujabc — ^a^iibc- Note the appearance of the spin con- 
nection, required by covariance, and the powers of g next to the momentum operator that 
ensure hermiticity of the hamiltonian. This hamiltonian is general enough to ahow future 
apphcations to the first-quantized theory of higher spin fields on curved backgrounds. 
The corresponding euclidean classical action in configuration space is given by 



where D^ip^ — dri^^ ^ x^oJiJ^bi^i- describes the particle propagation for a euclidean time 
/3 and will be used in the path integral quantization. For notational simplicity we do not 
make explicit distinction between quantum operators and classical variables, as it will be 
clear from the context which one is used. 

Having defined the model, we proceed with the rest of the paper and start describing 
the computation of the transition amplitude. In section 2 we use operatorial methods. 
In section 3 we consider path integrals in various regularization schemes, namely TS, MR 
and DR, calculating the corresponding counterterms and finding complete agreement with 
the expression of the transition amplitude found in section 2. Finally, we present our 
conclusions and perspectives in section 4. 

2 Transition amplitude from operatorial methods 

The aim of this section is the explicit computation of the transition amplitude for a eu- 
clidean time /3 determined by the quantum mechanical hamiltonian H given in (1.1). We 
compute the matrix elements of the evolution operator e~^^ between position eigenstates 
for the bosonic variables and suitable coherent states for the fermionic ones. We consider 
a short time expansion, and using systematically the fundamental (anti)-commutation re- 
lations of the basic operators x^p and '0 we obtain the final perturbative answer. The 
calculation is tailored after similar ones performed in [9] for the N — 0,1,2 supersym- 
metric quantum mechanics, explained carefully in [10] for the bosonic case, and in [28] for 
a class of supersymmetric sigma models on Kahler spaces which arise from the study of 
certain higher spin equations on complex manifolds [29]. As typical in semiclassical ap- 
proximations, the result can be cast as the product of three terms: i) the exponential of 
the classical action evaluated on the (perturbative) classical solution, ii) a standard leading 
prefactor depending on the propagation time as /3~^ , usually arranged in the so-called Van 
Vleck determinant, iii) a power series in positive powers of the propagation time, which 
identify the heat-kernel coefficients [8]. In our calculation we keep the approximation up 
to the first non trivial heat kernel coefficient, i.e. up to order /3 in the power series, keeping 
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in mind that the bosonic displacement = — can be considered of order as 
explained later on. 

In this section we restrict the calculation to even N = 25^, with S an integer, which 
allows us to introduce complex combinations of the fermionic operators and consider the 
corresponding coherent states as a (over)-complete basis of the associated Hilbert space. 
This is appropriate for applications to fields of integer spin. The method can be easily 
extended to the case of odd N ^ appropriate for applications to fields with half-integer spin: 
one way is to introduce another set of auxiliary Majorana fermions ("doubling trick") to be 
able to consider complex (Dirac) fermions and their coherent states. However, we refrain 
from doing that at this stage, as we employ it in the path integral approach of the following 
section. 

Thus we consider N — 2S, and introduce S Dirac worldline fermions out of the 2S 
Majorana fermions i/jf 

= + 7 = 1,..., 5 (2.1) 

SO that the new variables display the usual creation-annihilation algebra 

With these complex fields at hand we can readily construct coherent states such that 
^/l^) — Vi\v) (A|^^^^ = (A|A^^, with usual normalization (AI77) = e^'^ (other properties 
of fermionic coherent states can be found in appendix A). Denoting by \yr]) = \y) ® \r]) 
where \y) is the usual position eigenstate, the euclidean heat kernel we want to compute is 

(xA|e-^^|yr?) (2.2) 

where the quantum hamiltonian (1.1) and the covariant momentum, written in terms of 
the new complex variables, read 

with = In close analogy with the procedure employed in [9] we first focus 

on the mixed amplitude containing momentum eigenstates on the right hand side 

(xAle-Z'^br;) = 5]i^/3'=(xA|i7'=br?) . (2.4) 

k=0 

From the above expansion in /3 one could naively think to retain only the linear term, 
e~^^ = 1 — /3H + 0{/3'^)^ to get the answer up to the desired order. This is not the case 
for nonlinear sigma models, as well-known (see [9]). In fact it is necessary to take into 
account contributions for every order k, the only truncation being the number of (anti)- 
commutators to keep track of. In order to evaluate (2.4), we push all p's and ^^'s in each 
factor to the right hand side and all x's and ^^'s to the left hand side, taking into 
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account all (anti)-commutators, and substitute with the corresponding eigenvalues. Since 
the hamiltonian is quadratic in momenta, the matrix element of is a polinomial of 
degree 2k in p, so that in general one finds structures of the form 



2k 



(xA|i?^|p77) = 5^5f(x,77,Ay(x|2.)(A|77) , (2.5) 

where stands for the part of such polinomial with precisely l-th powers of p eigenval- 
ues, and the coefficients Bf^ are computed in appendix B. For plane waves we use the 
normalization 

{x\p) = {27v)-^/^g-^/\x) e^^-^ 

so, inserting in (2.2) a momentum completeness relation, 1 = / d^p\p){p\^ and rescaling 
momenta by = g^, we obtain for the transition amplitude 



e 



oo 



(^k 2fc (2.6) 

k=0 ' 1=0 

Let us look at this formula: it is well-known that the leading free particle contribution to 
(2.2) contains the exponential exp{— }, that gives the effective order of magnitude 
{x — y) ^ V^- Therefore, we see that q ^ and it is immediate to reahze that, for all 
k, only the terms B^j^, ^2k-i -^2k-2 ^^^^ contribute up to order /3. Taking for Bj^ the 
expressions given in appendix B, it is immediate to sum the series in k 



+ /3 [ - ^9^. Ingd^ \ng- ^d^d^ \ng 



- ^g^d^ In 5 - {^d^^g^ + ^g^g^ + ^g^g^ + ^^r") <?m Qu + (^5'^""' + l^^^^g" 
+ Y^f'^f: + YS^r'') lo qx q, qu - {^^''^f''') <?A q, qu q. q. qp 

- \d''{g^'ux)q. qa - Ig^^^oo^ qxqa - Ig'^'oo^ f q.qa + ^5^"^^/. 5^^" quq^qaq^ 

- ^g'^'Lo^d, \ng+ ^d^i^g'^'^co,] + ^{g'^'cv^abCOucd - 8a RabcdK)^"- rfrf^ 

+ A«- ry^'A'^- ry'^) - \g^'' u^^^hg^^ u^.ai^''- v"^'' + X"- r?''A<^- v'')q.q, -v]] , (2.7) 

where we remind that all the geometric quantities are evaluated at the final point x, unless 
otherwise specified, and we used the following compact notations 



Now we are ready to perform the q integration, that reduces to a bunch of gaussian integrals. 
Defining the coordinate displacement as = y^ — x^, one obtains the transition amplitude. 
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expanded up to order /3, in a form that hides manifest covariance 

(xAle-'^^lyr?) = (27r/3)-^/2 [g{x)/g{y)\ 1/4^-^9^.-''-'' gA^.r," ^ g-i/i q^^i/a 



^ dxg^,. z^z-z^ + \ z^z^ g-^'^ d^d^g^l^ - ^ z^g-^l^d^g'l^dxg,, z^z^z 



'n9 
p 



1 

+ 2 



12/3 



+ ^ i?^, z'^z'^ + zi'iOf, + ^z'^z'S^cu, + ^zi'oo^ z'g^'^d^gx. 



V4/3 



5Ac 



+ /5 



(2.8) 



This form is quite exphcit. However, it can be simphfied and cast in alternative and more 
suggestive forms. Keeping in mind that the exponent of the on-sheh action should appear 
in the final result, one can factorize (2.8) in the following way^ 



(xA|e-'^^|yr?) = (27r/3)- 



-D/2 



exp {- z^z- + \ d^g,x z^z^z^ + ^ (d^d.gx. - \ g^r V% z^z-z^z"\ } 



exp {a'^- T]a + Z^^^a}X- lf + \z^ z"" (d^UJ.ab + CU^acCu/fc) A'^- rf - RabcdX"- V^X"' V"^} 



1 



12 



1 



12 



(2.9) 



The first exponential is the expansion up to order /3 of the on-shell bosonic action evaluated 
with the boundary conditions x^(0) = and x^(/3) = 



Sx 



f 

Jo 



/3 



dr 



on shell 



- g^, z^z- + - d^g^x z^z^z^ + ( d^d^gxa - 7; gpr T^, F^, ) z^z^z^z 



yll yl^ y\ yCF 



(2.10) 



where all functions in the second line are evaluated at the point x. Similarly, one can see 
that the second exponential in (2.9) is the expansion, up to order /3, of the fermionic action 
evaluated on-shell, with boundary conditions ^^j(0) = v(\ and ^^^^(/3) = A^^, and with the 
usual boundary term added 



dT 



on shell 



A«- r?„ + z^cj^afeA'^- 77" + -z^'z" ( di^uj.ab + oJi^ac^.'^b ) A"- v" " 4a/3 i?«6,dA«- 17" A^- r?' 



b \c ^ 

(2.11) 



^The factor [g{x) / gd/)]^^"^ cancel against its inverse, whose Taylor expansion around x can be factored 
out from (2.8). 
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where the covariant time derivative reads D^^^j = + x^cj^^^^^j. Similar calculations 
of on-shell actions can be found, for instance, in [9]. Once the expansions of the on-shell 
classical actions have been recognized, the transition amplitude can be cast in the following 
covariant form 



{x\\e-P"\yi^) 



1 



(27r/3)^/2 
1 



exp{-(5a; + 5^)} 
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12 



(2.12) 



with all functions evaluated at point x. This is our final result for the transition amplitude. 

For comparison purposes, it may be useful to present the result after tracing over the 
fermionic Hilbert space 



{x\TT^ie-^^)\y) 
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exp{-5a;} 



1 + —R^^uzf^z" + f3 



(27r/3)^/2 
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12 
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N 



16 

l + 6aiV 



.yH^v, , , ah 



(2.13) 
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R-V] +0(/32 



Evaluated at coinciding points {z'^ — 0) it reads as 



(x|T¥^(e-^^) 



ND 

2— 



(27r/3)^/2 



1 + /3 



l+_QaN 
12" 



R-V) +0(/3" 



(2.14) 



which shows the first heat kernel coefficient at coinciding points for our model. As we shall 
see, the last two formulas remain valid also for odd N . 



3 Transition amplitude from path integrals 

In the present section we compute the transition amplitude by making use of path integral 
methods. To define the path integrals we fix a suitable regularization scheme, and compute 
the transition amplitude. Then, by comparing with the previous section or alternatively 
by imposing the Schrodinger equation, we identify the corresponding counterterms that 
ensure scheme independence. 

Unlike the previous section here we treat together both the cases with even and odd 
numbers of Majorana variables. In order to do so, we found it convenient to use the 
so-called "doubling trick" that consists in doubling the number of fermionic variables by 
adding "spectator" Majorana fermions il^'^ that satisfy free anticommutation relations. 
These new fermions are spectators in that they do not enter the interactions. With the 
help of these new variables one can consider Dirac fermions 

*i =^(V'f + #r), ^'? = ^(v?-#f) (3.1) 

that satisfy 

{*f,*5} = r7%- (3.2) 
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giving rise to a set of fermionic harmonic oscillators, whose Hilbert space can be represented 
in terms of bra and ket coherent states 



^t\v)-vt\v), (A|#f = (A|Af (3.3) 

and whose properties are briefly recalled in appendix A. The wave function of the system 
$(x,A) = (a;A|$), with {xX\ = {x\ ® (A|, evolves in euclidean time as 

$(x,A;/3) = {xX\e-l^"\^) 

= J d^'y^M J dCdv e-f-^(xA|e-'^^|yr?)$(y,C;0) (3.4) 

where we have used the spectral decomposition of unity for the bosonic and fermionic 
sectors 

l6 = y d''y^/9{y)\y){y\ , if ^ J dQdrj e-^"^|r?)(C| . (3.5) 

The evolution is generated by the basic transition amplitude that can be written in terms 
of a path integral as 

(xA|e-^^|y77) = (xA;/3|y7/;0) 

= / DxDaDbDc / D^D^ ^-S[x,aAc,^M ^3 

Jx(-l)=y J^{-l)=J~Bri 



where 



S[x, a, b, c, ^,^] = ^ t dTg^Mr)){x''x^ + a^'d' + h^d') 
J-1 



+ ^ I dTRabcd{x{T))i,'ii,\i,'li,'^ + ^ I dTV{x{T)). (3.7) 



Here, and in the following, we use a shifted and rescaled euclidean time r G [—1,0] to 
make comparison with the literature easier. We use bosonic (a^) and fermionic (6^, c^) 
ghosts to reproduce the reparametrization invariant measure Vx — Y\r V9i^i^))^^^(^)- 
We also rescaled fermionic "coordinates" and ghosts so that all propagators are of order 
/?, and added a fermionic boundary term to be able to set boundary conditions at initial 
time for and at final time for ^. Finally, let us note that the arbitrary potential V will 
eventually be modified by the addition of a counterterm Vct related to the regularization 
chosen. Apart form these modifications, and the inclusion of the spectator fermions, this 
is the same action given in eq. (1.3). 

We are interested in the short-time perturbative expansion of the transition amplitude. 
Thus we expand all geometric expressions about the final point x^, and split the action 
into a quadratic part plus interactions, 5 = 5^2 + Sint^ with 

^2 = j^g^u dr {x^x^ + a^a- + V^d") + dr ^>aM - -^^aA^m ■ (3-8) 
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Here and in the following, geometric quantities with no explicit functional dependence are 
meant to be evaluated at the final point x(0) = x, such as ^^^y = ^^^^(x). We can thus 
split the fields into backgrounds, satisfying the free equations of motion with corresponding 
boundary conditions, and quantum fluctuations; namely 

x^{r) = x^(t) + q^{r) , x^{r) = - z^r , q^{0) = q^{-l) = (3.9) 

^Kr) = ^Kr) + QKr) , ^^r) = ^/Pvt Q^-^) = (3.10) 

^^(r) = ^{r) + Q^(t) , ^(r) = , 0,^(0) = (3.11) 

where ^ — . The free on-shell classical action reads (henceforth we use a dot to 
indicate the contraction of whatever type of free flat indices) 

S2 = j^9f..z'^z''-X-7j (3.12) 

and the free propagators for the quantum fields, together with their Feynman diagrams, 
are 

(g^(r)g'^(cT)) = -f3g'''' AiT,a) = (3.13) 

{a^^{T)a''{a))=l3g^''Agh{T,a)= (3.14) 

{b^{ry{a)) = -2pgf^''Agh{r,a) = (3.15) 

{QUr)Q'j{a)) = /3ry%G(r,c7) (3.16) 

where the right-hand sides are given, at the unregulated level, by the following distributions 

A(t, ct) = T(a + 1)9{t - ct) + a(T + l)9{a - r) 
Agh{T,a) = S{T-a) (3.17) 
G(r,a) = -9(a-r) 

which obey the Green equations **A(t, a) = Agh{^i cr) = ^(^ — ^) and *G(t, a) = 5{t — a), 
with boundary conditions A(0,cr) = A(t,0) = A(-l,a) = A(t, -1) = and G(0,a) = 
G(t, —1) =0. Dots to the left (right) indicate derivatives with respect to the first (second) 
variable. These propagators have to be regulated in order to deal with divergences and 
ambiguities present in some diagrams. However, for each regularization scheme, one is able 
to cast all expressions in a way that can be unambiguously computed by using directly the 
expressions (3.17). 

Since fermions enter the interactions only through the combination = "73 '^^^ ^ 
it is convenient to write backgrounds and propagators for these fields as well. We split 
them as '0f (r) = '^/Sf + x^{^) with 

^? = /f(^r + A?) (3.18) 

and 

(X?(r)x5(c7)) =/3r7'^^5i,Af(r,<7)= (3.19) 
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satisfying •Ai?(T, a) = — A^(t, a) = 6{r — a) and given at the unregulated level by 



A^(r, a) = l{e{r -a)- e{a - r)) = ^6(r - a) 



(3.20) 



We only wrote propagators for unprimed fermionic fields since only they enter the interac- 
tions. The primed fields instead only contribute to the one-loop semiclassical factor that 
normalizes the path integral and drop immediately out of the computation. 

The interaction vertices that enter the perturbative expansion can be obtained by 
Taylor expanding the action (3.7) about the final point x(0) = x and read Sint — Yl^=z 
with 



(3.21) 



(3.22) 



+^Rabcd dT{r + x") • {v + xiir + x") ■ ir + Xl 

Higher order terms are not reported because we compute the transition amplitude to order 
/3, for which only the previous two terms are needed; all fields, classical and quantum, 
count as 

The transition amplitude can now be computed perturbatively using Wick-contractions 
of the quantum fields 



(3.23) 



with the suffix c indicating connected diagrams only and with the prefactor A soon to be 
commented upon. The various Wick-contractions produce 



exp 



1 



(53)-(54) + -(53^) 



exp 



^ dxg^.^z^z^z" + ■ 



4/? 



2/3 



'xOag^^.z'^z'^z^z^ - ^d^gl^ - ^d^d'^gxaU 



+ \dxdagz^z'' I5 + \d^gx.z^z'' le + dxgaz^z'' h + j^dxu^abr ■ i^'^z^z'^ 
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- ^{dxg,^ufl9 - ^dxg^ud^gxuho - ^{dxgflu - ^gxdxglu - ^gl^is 
+ ^dxg,j.udxg^'i,z^z^' lu + ^dxg^udi,g^^xz^z^' lis + ^dxg^udyg^uZ^z^' lie 
+ dxg^jiud^g^'^z^z^' Ii7 + ^dxg^udxgz^z"" lis + ^dxg^ugxz^z"" I19 

+ ^dxg,^iyd,^gz^z''l2o + dxg,^i,g,^z^z''l2i - ^dxg^udxg^'u'Z^z'' z^' z""' I22 

- ^dxg^vdwgxu'z^z^'z^' z""' I23 - ^Oxg^vdyg^j^'Z^z'^z^' z""' I24 

+ ^^Aa6l25 - ^UXab^a^'^Z^z'' I26 + ^C^Aac ^A6^< ' l27 

- ^CJAac^a^"^^^^^^ • ^^l28 - ^^Aa6^Acd^^ ' ' I29 

+ ^^Aafe^A^'' • I30 + ^^Aa6^A^^ ' I3I " ^^Xahdxg^^uZ^ Z^ i^"" • I32 

-^^Ma65A^M-^^^>''-^^l33 



(3.24) 



where we made use of several shortcut notations, including g^ — g^^dxg^u^ dxg = g^^dxg^u^ 
dxdag = g^^dxdag^u and = g^^dxda- The integrals 1^^ are reported in appendix C along 
with the pictorial representation of the diagrams they belong to (an integral named Is 
is absent, but we kept the same notation as in [10], where such an integral arose from 
the coupling to external gauge fields, to allow easy comparison). We compute them in 
the following subsections, using the different regularization schemes discussed in the intro- 
duction. The purely bosonic contributions (A; < 24) are well-known from many previous 
computations, see [10] for example; the remaining ones have been computed, though only 
for < 2, in the time slicing regularization technique [16], in dimensional regularization 
[25, 26], and in mode regularization [21]. However, in the last two schemes fermions were 
traced out to obtain directly heat kernel coefficients and trace anomalies. In the present 
case we are interested in the full transition amplitude and wish to keep N generic, so that 
we need to reconsider all such integrals with fermionic contributions. Finally, the prefactor 
A can be fixed by requiring that, in the limit /? ^ 0, the transition amplitude reduces to 
{x A; f3\y rj; 0) 6{x — y) e^'^, which in MR and DR gives 

In TS such prefactor can be deduced directly starting from operatorial methods, and reads 



A 



9{y) 



1/4 -^ 

(3.26) 



(27r/3)-D/2 • 

We are now ready to consider the various regularization schemes. 
3.1 Time slicing regularization 

Time sliced path integrals in curved space were extensively discussed in [15, 16]. In essence 
time slicing regularization consists in studying the discretized version of the path integral. 
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as derived form the operatorial picture by using Weyl ordering and midpoint rule, and 
recognizing the action with the correct counterterms together with the rules how to compute 
Feynman graphs that must be used in the continuum limit. These rules state that the Dirac 
delta functions should always be implemented as if they were Kronecker delta's, using the 
value 0(0) = ^ for the discontinuous step function. We do not need to repeat here many 
computations, namely the bosonic ones, though one can easily compute them using the 
explicit expressions collected in appendix C, or extract them directly from [10]. It is 
enough to focus on the graphs containing fermionic lines, i.e. 1^ with k > 25. 

Thus, let us start with I25, that is the only diagram that depends heavily on the 
regularization chosen. Following the TS prescriptions, we have 

rO rO rO rO 

^25 



/ /* f f / \ / 1 \ 2 

/ drda •A^Al^ / / drda^l - 5{t - a)j (^-e{T - a)j 

= i drda e\r - a) = - j _^ drda = - (3.27) 

since 'A* = 1 — 5{t — a) and e(0) = 0, as follows form 0(0) = \ and eq. (3.20). The regular 
I26 does not need any prescription and corresponds to the last line above, giving I2Q — \. 
Similarly, one finds that the integrals 1^^ with k > 27 yield a vanishing result (those with 
k > 28 actually contain only bosonic propagators, but depend on the external fermionic 
backgrounds) . 

The transition amplitude then reads 

{x A; /3\y r/; 0) = (27r/3)^/2 (bosonic) + —uxabZ^i^"" • 

+ (3.28) 

where (bosonic)^^ contains the purely bosonic contributions of (3.24), including the metric 
factors appearing in (3.26). It can be extracted from the literature, or easily recomputed 
with the set-up described above, and reads 

(bosonic)^, = - ^ d^g.x z'^z'^z^ - ^ {d^d.gxa - \ Qpr V% n,) z^z^z^z'' 

+ ^^R,.z^z^ + /3(^5^^r^,r^, _ _Li? _ . (3.29) 

This transition amplitude is the one computed with a TS regularization of the Feynman 
diagrams, and in general differs form those computed with other regularizations if no 
counterterms are introduced. In particular, eq. (3.28) satisfies a Schrodinger equation 
with a non-covariant hamiltonian that differs from the one given in eq. (1.1). One expects 
different regularizations to be related by finite local counterterms, so we need to identify 
the correct counterterm that make sure that we are discussing the same quantum theory as 
the one associated to the hamiltonian i7 of eq. (1.1). To achieve this we can either compare 
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with the transition amphtude obtained by operatorial methods, or compute directly the 
hamiltonian associated with the transition amphtude above. We shah fohow both methods 
as a check of our calculations. 

As it stands the transition amplitude calculated above cannot be compared directly 
with the result from canonical methods in eq. (2.12): the latter is valid for even N and 
is written in terms of fermionic coherent states that correspond to the physical fermions 
only, while the present result, valid for any integer N, contains also the Majorana spectator 
fields introduced in the fermion doubling trick. To overcome these differences, we can take 
a trace over the fermionic Hilbert space and eliminate the contribution of the decoupled 
spectator variables by subtracting their degrees of freedom. Thus, let us compute the trace 
in the fermionic sector of eq. (3.28) by integrating over the Grassmann variables with 
measure J dr]dXe^'^^ see appendix A. The final result can be obtained, for instance, using 
the standard Wick-contractions associated to the gaussian integral J drfdX e^^'^ (note the 
factor 2 in the exponent arising form the trace measure and the leading part of (3.28)), 
which gives the following propagators 

(A«7?5) = ^ry% ^ W^'j) = (3.30) 

where we used the definition in (3.18). Note that the normalization factor J drfdXe^^'^ = 
has to be renormalized to 2^ to undo the fermion doubling. Taking all this into 
account, and setting z = — x = for simplicity, we obtain 



(3.31) 

where the first term in round brackets is due to the purely bosonic contributions. Com- 
paring with (2.14) one recognizes the counterterm that needs to be added to the potential 
V in the path integral action to achieve equality 

^^s' = -{1 + ^)r + yrurl + ^^.^uT- . (3.32) 

Power counting shows that no higher order corrections to the counterterm are to be ex- 
pected. 

Alternatively, one may compute the hamiltonian appearing in the Schrodinger equation 
satisfied by the amplitude (3.28). To do this we insert the latter into (3.4), Taylor-expand 
to order /? all terms in the right hand side about the final point, and identify the Schrodinger 
equation associated to it. Comparing with the Schrodinger equation due to the hamiltonian 
(1.1) one deduces eventual counterterms. Let us describe this computation. We perform 
the Gaussian integrals over d^y = d^z and the integrals over the fermionic coherent states 
using the properties summarized in appendix A. The purely bosonic contributions of the 
diagrammatic expansion yield the standard TS result 

$(x. A; /?) = f 1 - /3dt - /3Hb + 0{l3^))^x, A; /3) (3.33) 
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with 



Hb = -^d^g^'^Vgdu + V + -^R- g'^'T'^^x^t,) . (3.34) 

that by itself requires the addition of the famihar counterterm Vr^^^ — —^(^R — g^^T^^T^^^ 

into the path integral, in order to get the covariant H = ~2^^m^^^V^^^ ^ need 
to dress this result with the fermionic contributions. One way to perform the integrals 
over the fermionic coherent states is to use the formulas involving fermionic delta functions 
derived at the end of section A. By defining the full hamiltonian Hp as 



Hf = Hb + A (3.35) 



and using (3.33), we obtain 



+aRabcdM'''M'^^ + - ^{coxab? (3.36) 

z lb 

with M^^ = ^ f A'' • + A'' • ^ - • ^ + ^ • ^] being Lorentz generators. We observe 
that the noncovariant terms in the line above are those necessary to render the bosonic 
hamiltonian Hb local-Lorentz covariant, namely 



1 

where i?^ is the fully covariant derivative (with both Lorentz and Christoffel connections), 
so that in order to have^ 

H = -^g^^-D^D, + aRabcdV ■ V'V' ■i^'^ + V (3.38) 

one needs to add to the path integral the same counterterm found before in eq. (3.32). 
Thus we found complete agreement. 

The above expression for the counterterm V^^^ matches all the previously known 
cases [15, 16]: the purely bosonic case {N = 0) is obviously reproduced. For = 1 
supersymmetry fixes a = —1/4 and y = so that V^^^ = ^g^^F^^F^^ + ^(^/iafe)^ comes 
out correctly. Note that in this case the relation Rabcd^^^^^^i^^ = allows to use 

a = and V = |i? as well to identify the same supersymmetric hamitonian. For N — 2^ 
supersymmetry requires a = — | and y = 0, so that again V^^ = ^^^^F^^F^^ + \ {oo,jiahY 
is correctly reproduced. 



Hf = Hb + A = -^gf^-D^D, + aRabcdM'^^M'"^ + V 



^Here we use (A|V'"*|$) = (^gfr + $(A), so that the Lorentz generators can be written as M°-^ = 
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3.2 Mode regularization 

In this section we approach the previous computation using mode regularization (MR), 
that can be considered as the worldhne equivalent of a cut-off regularization in momentum 
space of standard quantum field theories. Mode regularization starts by expanding all fields 
in Fourier sums, thus identifying a suitable functional space to integrate over in the path 
integral. The regularization is achieved by truncating the Fourier sums at a fixed mode 
M, so that all distributions appearing in Feynman graphs become well-behaved functions. 
Eventually one takes the limit M ^ oc to obtain a unique (finite) result. Often one can 
proceed faster: performing manipulations that are guaranteed to be valid at the regulated 
level one may cast the integrands in alternative forms that can be computed directly in 
the continuum limit, without encountering any ambiguity. 

The boundary conditions for the bosonic variables are as in (3.9), so that the bosonic 
quantum fluctuations, as well as the ghosts, are naturally expanded in a Fourier sine series 



where the mode M is the regulating cut-off that is eventually sent to infinity, as in [17]. 
This choice preserves the boundary conditions imposed at initial and final times. On the 
other hand, fermionic fields satisfy first order differential equation and carry boundary 
conditions only at initial or final times, but not at both, see eqs. (3.10) and (3.11). Thus 
we find it useful to expand the fermionic quantum fields in half integers modes as follows 



This choice preserves the boundary conditions and provides us with a regulated functional 
space to integrate over also for the fermions. In addition, it produces a simple regulated 
kinetic term that is easily inverted to obtain the propagators. Finally, the path integral is 
defined as a regulated integration over the Fourier coefficients of the various fields. 

Perturbatively, the propagators are as in eqs. (3.13)-(3.16) and are regulated as follows 



M 




(3.39) 






m=l 



M 




(3.42) 



m=l 





As a consequence 




(3.44) 
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for (3.19) that turns out to be translational invariant. Note that in the hmit M oc the 
previous formulas reproduce eqs. (3.17) and (3.20). 

We are now ready to compute the Feynman integrals with fermionic contributions in 
MR (the purely bosonic ones are standard and can be extracted from [10]). Using the 
regularized expressions one obtains I27 = I28 = ••• = I33 — 0. Also, one finds I26 = |, 
as it is unambiguous and gives the same result in all regularization schemes. The only 
tricky integral is I25. Using the mode regulated propagators we can manipulate it using 
an integration by parts as follows 



^25 



/O rO 
-1 J-1 



drda •A^Al 



dr 



a=0 



rO 



drda A** 



dr 



A'Al 



1(7=0 



1 

1 + 2 



rO 



(3.45) 



drda CA* - A**) A^ 



The first term is unambiguous and can be computed directly upon removing the cut-off to 
give |. In the second term one can use the following identity valid at the regulated level 
(with X — T — a) 



•A-(r,c7)-A"(r,c7) 



-2cos(-) 



•Ai7(x) + 1 - cos(7rMx) . 



(3.46) 



The integral involving the last term (cos(7rMx)) produces a vanishing result in the limit 
M ^ oc, thanks to the Riemann-Lebesgue lemma, while the rest can be computed using 
an integration by parts and removing the regularization to yield the final result 



l25 



drda 



cos 



(t) 



1 



(3.47) 



The transition amplitude computed in mode regularization hence reads 

1 



{x\;/3\yri;0) 



(27r/3)-D/2 



■ exp 



(bosonic) 



I3N 2 



(3.48) 



where (bosonic) J^^^ is the purely bosonic contribution that can be found in [10]. Comparing 
the above result with the TS one given in eq. (3.28), we notice that the only difference 
coming from the fermionic sector sits in the coefficients of the last uP' term, which is due 
to I25. As the bosonic part of the MR calculation requires the counterterm 



24 \ v\) 



to reproduce the heat kernel for H = -j^d/^g^^'^du, where (r^;^)^ = g""' g^^' g^^/T'^^T'^'y, 
we see that, in order to obtain the correct amplitude for (1.1), we need the counterterm 



(N) 



MR 



aN 



N 



24 



(3.49) 



that again matches the results known for A'' = 0, 1, 2 [19, 21]. 
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3.3 Dimensional regularization 

Finally, we reconsider the previous calculations using dimensional regularization. This is a 
perturbative regularization that uses an adaptation of standard dimensional regularization 
to regulate the distributions defined on the one dimensional compact space / = [0,/3]. 
One adds d extra non-compact dimensions and perform all computations of ambiguous 
Feynman graphs in d + 1 dimensions. Extra dimensions act as a regulator when d is 
extended analytically to the complex plane, as in the usual QFT case. After evaluation 
of the various integrals one should take the d — > limit [24]. This is quite difficult in 
general, since the compact space / produces sums over discrete momenta, and standard 
formulas of dimensional regularization do not include that situation. However one may 
use manipulations valid at the regulated level, like differential equations satisfied by the 
Green functions and partial integration, to cast the integrands in equivalent forms that can 
be unambiguously computed in the d ^ limit. While purely perturbative, this method 
carries a covariant counterterm that simplify extensive calculations, as the one performed 
in [30] to obtain trace anomalies for a scalar field in six dimensions using worldline methods. 

In DR the computation turns out to be quite simple for most diagrams: I27 = I28 = 
because the integrand is odd, whereas I29 = • • • = I33 = as can be shown by integrating by 
parts. Also I26 = |, as it is regular and can be safely evaluated by using the unregulated 
expression for the propagator. As usual, more care is needed to compute I25 since the 
integral is ambiguous (products of distributions). By dimensionally extending the cubic 
vertex 



where A{t,t') = -^'^ {q{t)q{t')) and AF{t,t') = /3"^(x(^)x(^0) (the bar indicates Ma- 
jorana conjugation, see [24]) are the dimensionally regulated propagators where — 



-Ap^'^ 5(^+^)(t,t'), and AAsit^f) = dA&BA{t,t'). One can partially integrate the above 
integral without picking boundary terms since in the compact direction bosonic fields van- 
ish at the boundary, and in the extended directions Poincare invariance allows partial 
integration as usual. Hence the above expression becomes 



(3.50) 



the above integral becomes 




(3.51) 




j dtdt' Ab 9^tr[7^AF7-^AF] 
j j dtdt' Ab tr {^Af)i^ Ap - Api^iAp^') 
2 [ dt Ab tr[-f^AF{t,t)] 125 = 2/" dr 'A A^lr = 



(3.52) 



-1 



where we have used the regulated Green equation, integrated the delta function, and used 
that Ai? vanishes for coinciding points. 



The transition amplitude in DR thus reads 




(3.53) 



where the bosonic contribution (bosonic)^^ can be extracted from [10]. Comparing this 
result with the ones obtained before in other regularizations we note that, thanks to the 
vanishing of I25, the diagrammatic expansion does not produce the term proportional to 
cj^ which previously had to be canceled by the counterterms. Thus, the standard bosonic 
counterterm V^l = — ii? is dressed up to 



This matches known results for = 0, 1, 2 [24-26]. 
4 Conclusions 

In this paper we have studied the quantum mechanics of one dimensional nonlinear sigma 
models possessing a O(A^) extended supersymmetry on suitable target space backgrounds. 
Considering an arbitrary background geometry, we have computed the transition amplitude 
for short propagation times using both canonical and path integrals methods, obtaining in 
the latter case the correct counterterms associated to various regularization schemes which 
are needed to define unambiguously the path integrals. 

A possible use of our results is in the discussion of higher spin fields in a first quantized 
picture. Worldline approaches are useful in finding efficient ways of computing amplitudes 
for relativistic processes both in flat space [31] and curved spaces [32]. The quantum 
mechanical nonlinear sigma models discussed here arose precisely in an attempt to use 
worldlines methods to describe one-loop effective action due to higher spin fields in a curved 
background [7, 27, 33]. In future works we plan indeed to use the path integrals constructed 
here to study effective actions induced by higher spin fields and compute corresponding 
heat kernel coefficients. Finally, one might wish to extend our results to the OSp(A^, 2M) 
quantum mechanics of ref. [34], which have also found applications to higher spin theories 




(3.54) 



[35-37]. 
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A Fermionic coherent states 



The even-dimensional Clifford algebra 

{^^, = 6^^ , M, TV = 1, . . . , 2/ (A.l) 

can be written as a set of / fermionic harmonic oscillators (the index M may collectively 
denote a set of indices that may involve internal indices as well as a space-time index), by 
simply taking complex combinations of the previous operators 

= + (A.2) 

al = ^(^^^-i^^^') , m = l,...,; (A.3) 
{a-,a],} = C (A.4) 

and it can be thus represented in the vector space spanned by the 2^ orthonormal states 
l'^) — nm(^J^)^'^lO) with ajn\0) — and the vector k has elements taking only two possible 
values, km = 0, 1. This basis (often called spin-basis) yields a standard representation of 
the Clifford algebra, i.e. of the Dirac gamma matrices. 

An alternative overcomplete basis is given by the coherent states that are eigenstates 
of creation or annihilation operators 

(f^l = (0|e^-«'" ^ {f]\al = (f?|f?^ = fimifil ■ (A.6) 

Below we list some of the useful properties satisfied by these states. Using the Baker- 
Campbell-Hausdorff formula e^e^ — e^e^e'^'^l, valid if [X, y] = c-number, one finds 



(r?|0 = e^-« (A.7) 



that in turn implies 



(f?|a-|0=r(r?IO = (A.8) 

OTjjYi 

{v\al\C)^f,mm (A.9) 
so that {-^,fin} ^ 5'^. Defining 

df] ^ dfji ■ ■ ■ df]i , ^ d^^ ■ ■ ■ d^'' (A.IO) 
so that df]d^ = df]id^^ dfj2d^^ ■ ■ ■ dfiidS} , one finds the following relations 

dfid^ e-*'-^ = 1 (A.ll) 

dfidi e-^-^ = 1 (A.12) 
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where 1 is the identity in the Fock space. One can also define a fermionic delta function 
with respect to the measure (A. 10) by 



S{f] -X) = - A^) . • • {ff - ^0 = y e(^-^>^ . (A.13) 
Finally, the trace of an arbitrary operator can be written as 

Tr^ = y « e-^-«(-^|A|0 = j didf) e'^-«(f?|^|0 • (A.14) 
As a check one may compute the trace of the identity 

Trl = didf) e^'^{f)\i) = j d^df] e^^< = 2^ . (A.15) 

Let us end this section by deriving a few expressions involving fermionic delta functions 

that are helpful in the computation of section 3.1. Using ijj^ = \f^iVi + K)^ defined in 
(3.18) (here fermions are labelled by two indices, a tangent space index a and an SO(A^) 
R-symmetry index z), we may compute 



/ 



+ K) dCSiC-p) ^ (A.16) 



P=\ 



where the above apparently baroque notation is used because the fermionic derivatives only 
act upon the delta-function and not on eventual A-dependent expression that may appear 
next to it. We thus have 

J dCdr] e^^-f^V" • = f^M"'' J dC 6(C - p) _^ (A.17) 



with 



M"''^- A^-A^ + A". — -A"- — + — .— . (A.18) 
2V dpb a Pa dpa dpb/ 

One can then switch to the Lorentz generators 

by suitably subtracting those terms that appear when derivatives on A act on eventual 
functions of A that may show up before the wave function, such as M^^ in (3.36). 

B Bi coefficients 

The coefficients Bj^{x^f],^)^ defined in (2.5), can be computed following the strategy de- 
scribed in detail in [9] for = 0, 1,2, and in [28] for the complex U{N\M) sigma model. 
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First of all we divide the hamiltonian (2.3) into three pieces contributing at most two, one 
or no p eigenvalues 

Hb + Hi + H2 , where 
1 



H 

Hb 
Hi 

H2 



^9 -'i'vd''rv.9-'i' 



(B.l) 



First of all, notice that Hb is precisely the usual bosonic quantum hamiltonian, carefully 
studied in the literature [9, 10]. Let us start with B\^\ the only way to have 2k p eigenvalues 
is from k factors of Hb and no commutators taken into account, giving simply 



2k 



(B.2) 



For Bl^_ 1 we can have two kinds of terms. The first comes from k factors of Hb with 
one p acting as a derivative; this gives the corresponding of the purely bosonic 

model, whose computation is explained in detail in [9, 10]. The other term comes from 
k — 1 factors of Hb and one Hi, by substituting all operators with the corresponding 
eigenvalues. Putting things together we obtain 



TDk 2k-i 
^2k-lP 



ik 



k-i 



k-2 



1 



i k 



k-l 



(B.3) 

For B2j^_2 four types of term contribute: i) k factors of Hb^ giving the coefficient of the 
corresponding bosonic model, ii) k — 1 factors of Hb and one Hi, with one p acting as 
a derivative. This contribution gives four terms: the derivative acting from one Hb to 
Hi, from Hi to one Hb, within Hi or within the k — 1 Hb^^- Hi) k — 1 factors of Hb 
and one H2, substituting all operators with their eigenvalues, and iv) k — 2 factors of 
Hb and two Hi, substituting all with eigenvalues. Remember that in in) and iv) {^^,^} 
anticommutators have to be taken into account in order to obtain eigenvalues on the 
coherent states. Altogether it results in 

.2\ k-l 



TDk 2k-2 
^2k-2P 



1 



1 



1 



1 



—d^ \ng^^' \ng+ -S^S'^ Iny + -g^'d^ lug 



32 



k-2 



-^d^g^ + \g'g^ + -/g^"" + \gT\ 



Pl^Pu 



k-3 



k-4 



k-2 



rl 



+ {g^'^g" + Ig'^^'g^p" + -^g^d^"' 



Pi^PvP\Pa 



PuPvPWaPpPr 



d^[g 



1 



k-2 
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+ 



k — 1\ f p 



2\ ^-3 



9'''^aPp + \k(P^ 



2\ k-1 



9''''ujnd,\ng 



1 / 2\'=-l 



2 V 2 



+ 1- i^n"- i'^ 



2V\ 



2) (^) ^''"^'^"''^"''^A'^'^ 



PuPa (B.4) 

where again we use the compact notations for tensors introduced below eq. (2.7). 



C Feynman diagrams 

We hst the set of Feynman diagrams and the associated worldhne integrals that enter the 
computation of the transition amplitude to order /3 in section 3.. We are not reporting here 
those diagrams that involve fermionic self-contractions as with the rules of section 3 such 
self-contrations are trivially vanishing. Hence, the a priori non-trivial diagrams entering 
the contribution (53) are 



o — c 



(C.l) 
(C.2) 



Those contributing to {S/^) are 



00 ^OC> 



•A) I. 



dT 'A\ 




+ 




J\tt\'A' + "A)\, 





J dT A\r 

J dT rA'l 



(C.3) 
(C.4) 
(C.5) 
(C.6) 
(C.7) 



The remaining ones contributing to {S'^)c can be devided into purely bosonic contributions 
; + | I I dTda A{*A*^ (C.8) 



Iio = 



(C.9) 
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23 



and those with mixed bosonic-fermionic contributions 
I25 = ^ 



drda Al 



"-26 



1 J-1 

rO 

1 J-1 



drda 



"-27 



= <~> = 



y drda "A* Ai 



I28 



"-29 



I30 = 



L3I 



"-32 



"-33 



rO 



drda 



-1 J-1 
rO 



drda 



1 J-1 



rO 



drda 'A ('A* + "A)|^ 



1 J-1 



Q= drda 'A* A'l^ 



rO 



drdcr 'A 



-1 J-1 
rO 



drda a 'A* 



1 J-1 
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